Detectability of GRB Optical Afterglows with Gaia Satellite 
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ABSTRACT 

With the launch of the Gaia satehite, detection of many different types of transient sources 
win be possible, with one of them being optical afterglows of gamma-ray bursts (GRBs). Using 
the knowledge of the satellite's dynamics and properties of GRB optical afterglows, we performed 
a simulation in order to estimate an average GRB detection rate with Gaia. Here, we present 
the simulation results for two types of GRB optical afterglows, differing in the observer's line 
of sight compared with a GRB jet axis: regular (on-axis) and orphan afterglows. Results show 
that for on-axis GRBs, less than 10 detections in five years of foreseen Gaia operational time are 
expected. The orphan afterglow simulation results are more promising, giving a more optimistic 
number of several tens of detections in five years. 



Subject headings: Gamma-ray bursts 
1. Introduction 

European Space Agency's cornerstone mission 
Gaia is scheduled for launch in 2013. It is primar- 
ily an astrometric mission, designed to measure 
precise 3D positions for up to a billion stars in 
our Galaxy and other bright objects in the sky. 
In addition to astrometric measurements, spec- 
troscopy and multi-band photometry will be per- 
formed, providing radial velocities and additional 
astrophysical information on the stars and other 
objects. The main scientific goal of Gaia mission 
is to clarify the origin and formation history of our 
G alaxy. The goal s of th e mission a r e des cribed 
in lPerrvmanet all (l200ll).lLindegrenl (IpOlOt ). and, 
more thoroughly, in lTuron et al.l ()2005l ). 

Gaia will continuously scan the sky, covering 
the whole sky and building its own catalogue of 
celestial objects. A preliminary catalogue will be 
released every 6 months and a final catalogue, ex- 
ploiting all the data recovered during the mission, 
will be released at the end of the ope rational time 
(around 2020) (JLammers et al.ll2008l ). Comparing 
further observations with this catalogue will en- 
able Gaia to detect various transient sources, with 



some of them being cataclysmic variables, super- 
novae, active galactic nuclei, asteroids, etc. Since 
the catalogue of Gaia astrometric and photomet- 
ric data will be available only late into the Gaia 
operation, triggers from transient phenomena will 
actually be the first Gaia data released to the sci- 
entific community. In order to be prepared for this 
stream of data and to make sure the alert stream 
is accurate, reliable, and free from contamination, 
the Gaia Science Alerts Working Group has been 
formecQ. 

One type of possible transients to be detected 
by Gaia is gamma-ray burst optical afterglows. 
Gamma-ray bursts (GRBs) are extremely ener- 
getic e xplosions o c curring at cos mological dis- 
tances ( Piranll2005t iMeszarosI 120061 ). They repre- 
sent the most luminous events in the gamma-ray 
part of the electromagnetic spectrum known in the 
universe. Usually, the short-term prompt gamma- 
ray emission is followed by an afterglow emitted 
at longer wavelengths from X-rays to radio waves, 
which can last for several days. According to the 
standard fireball model (JRees fc Meszaroslll992 . 



-^ http://www.ast.cam.ac.uk/ioa/rcsearch/gsawg 



19941 ) the prompt gamma emission is produced 
in the internal shocks of the relativistically ex- 
panding ejecta, while the afterglow is synchrotron 
emission from relativistically accelerated electrons 
in the external shocks, created in the ejecta inter- 
action with the circumburst medium. In general, 
afterglow l ight curves exhib it the power-law decay 
with time (jSari et alJll998l ). althou gh many addi- 
tiona l features like density bumps ( Lazzati et al.l 
2OO2I: iGuidorzi et alj 120051). late energy inj ections 
(iNousek et al.l l2006t lEvans et all l2009l) and re- 



verse shocks (lAkerlof et al.lll999l:lFox et al.ll2003t 



Kobavashi fc Zhand I200T iGomboc et al.l I2OO8I ) 
have been seen in recent years, when better 
time coverage observations became possible due 
to accurate GRB a lerts from the Swift satellite 
(JGehrels et al.l 120041 ) and a number of follow-up 
robotic telescopes (e.g ., the Liverp ool Telescope 
(iGomboc et al.l l2005h. ROTSE (lAkerlof et al. 



20031 ). REM (jZerbi et al.ll200ll ). etc.). 

GRB explosions ar e believed to be collimated , 
rather than spherical (|Rhoadsl[l997l:lGranotll2007t) . 
An indirect argument supporting this is based on 
the high values of energy output in gamma rays 
(-Eiso), which is of the order of '^ Mqc^ if we as- 
sume that the explosion is spherically symmet- 
ric. It is known that GRBs occur in galaxies, 
and are of stellar origin (most probably the col- 
lapse of massive, rapidly rotating stars and merg- 
ers of compact objects). It is difhcult to conceive a 
model with a stellar progenitor, which can produce 
such high energies inferred in the case of spheri- 
cally symmetric explosions. More direct evidence 
comes from precise observations of afterglow light 
curves: an achromatic break in the light curve, 
which has been observed in many GRB af terglows 
at a f ew hours to days after the GRB itself (|Granot 
2OO7I and references therein), was predicted theo- 
retically prior to observational discovery (Rhoadg 



I999I : ISariet al.l[l999l ). The break (usually re- 



ferred to as a jet break) can be well explained 
by considering a jet with a half-opening angle 9j, 
moving initially at a relativistic speed with high 
Lorentz factor r '^ 100 — 1000. Because of the rela- 
tivistic beaming effect, the radiation emitted from 
the jet is beamed in a half-opening angle 9h = ^ 
in the direction of motion. The expanding jet is 
surrounded with a medium that causes the jet to 
decelerate. This results in an increasing value of 
the angle dh- Shortly after the prompt GRB emis- 



sion, we cannot distinguish between a spherical 
and collimated explosion due to strong beaming, 
since 9j > 9h- But after some time, the jet slows 
down and 9^, becomes larger than the opening an- 
gle of the jet itself. At this point, the radiation 
flux drops and we can see an achromatic break in 
the l ight curve (JKulkarni et al.l Il999t iKlose et al 



20041 ). As mentioned previously, a collimated ex- 



plosion results in a reduced estimation of energy 
output, i.e., Ej = ^9?Eiso, where Ej is colhmation- 
corrected energy in the case of a jet with half- 
opening angle 9y 

The shape of the afterglow light curve depends 
on the relative position of the observer's line of 
sight and the axis of the jet cone. We can de- 
scribe that with an angle ^obs between both di- 
rections. It is very unlikely for an observer to be 
looking precisely at the direction of the jet axis 
(to have 6'obs — 0). We get an 'on-axis' GRB and 
a 'regular' optical afterglow when 6'obs < ^j- If 
the jet cone is directed so that ^obs > ^j, an ob- 
server, due to the beaming effect, probably cannot 
observe the prompt gamma emission, i.e., detect a 
GRB. But the deceleration of the jet will gradually 
expand the beaming angle and the observer might 
detect the afterglow. The observed afterglow that 
was not preceded by a prompt gamma emission 
is called an orphan afterglow (OAJj. Even though 
OAs have not been conclusively detected yet, they 
have been studied a great deal over the last 15 
years (JGranotll2007l and references therein). 

As part of the scientific community preparation 
for the Gaia mission, we discuss Gaia^s potential 
for detecting these short, rapidly fading and un- 
predictable transients. We present here results of 
the simulation we performed in order to estimate 
the number of possible GRB optical afterglow de- 
tections with the Gaia satellite. The simulation is 
divided in two parts, the first one concerning on- 
axis GRB afterglows and the second one dealing 
with orphan afterglows. In Section 2 we briefiy 
discuss the dynamics and basic characteristics of 
Gaia. The simulation with all the necessary pa- 
rameters for both on-axis and orphan afterglows 
is described in Section 3. In Section 4, we present 
and discuss results of the simulation and possible 



^The OAs we are discussing here are not to be confused with 
optical afterglows, for which the observer does not detect 
the prompt gamma emission for some reason, even though 
the jot cone is directed into the observer's line of sight. 
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Fig. 1. — Sketch showing how Gaia wih scan the 
sky. The telescope axes are separated by an angle 
of /3 = 106.5°. The satellite spins with a constant 
angular velocity uj around its main axis, which is 
kept at an angle 7 = 45° from the Sun at all times. 
The main axis experiences slow precessional mo- 
tion around the Earth-to-Sun direction (17). Ro- 
tation around the Sun (angular velocity p) is also 
taken into account. Details of the dynamics are 
given in the Appendix. 



ways of afterglow identification. 

2. The Gaia Satellite 

For the purpose of our simulation, we need to 
know the way Gaia will scan the sky (JLindegren 
20101 ). Gaia will carry two identical telescopes, 
separated by the angle of /? = 106.5°, as shown 
in Figure [1] The satellite will make four rotations 
per day around its axis (which is perpendicular to 
the direction in which both telescopes are point- 
ing) with the constant angular velocity w. The 
direction of the axis itself is tilted by the angle 
7 — 45° from the direction of the Sun. The axis 
will experience slow precession motion around the 
Earth-to-Sun direction with a period of 63 days (il 
in Figure [1]). Gaia will have an orbit around the 
L2 point and will thus experience rotation around 
the Sun, which is shown in Figure [T] as a rotation 
around the x axis. Knowing w, fl, 7, /3 and p (one 
year orbital period around the Sun), we construct 



the scanning law of Gaia (see the Appendix R)) . 

Gaia's two telescopes will have a field of view 
of ^ 0.7° X 0.7° each. The expected hmiting 
magnitude in broadband G magnitude (details 
on the photomet ric system of Gaia are given in 
Jordi et al.lbOlOT) is G = 20 mag. The specifics 
of scientific performance are given on the official 
G aia World Wide W eb pagcl and, for example, 



m 



LindegrenI (|2010l ). There wih be a 7 



9 as- 

trometric CCD field in Gaia's focal plane. Each 
source will transit over the nine CCDs and will be 
observed by each of them with 4.4 s integration 
time. 

3. The simulation 

3.1. On- Axis Afterglows 

First we focus on on-axis afterglows, i.e., those 
with their jet cones turned in our line of sight. In 
such cases we can observe the GRB, which trig- 
gers a satellite and follow-up optical observations. 
Hence, we can base the initial parameters of our 
simulations on the actually observed GRB after- 
glow numbers and their characteris tics. Since the 



launc h of the Swift satellite in 2004 (jGehrels et al. 



20041 ) a large number of GRBs and their after - 
glows has been detected (jRoaming et al.l 120091 ). 
The Swift detection rate is about 100 GRBs per 
year. In about half of the detected GRBs there 
is no bright optical afterglow detected. Since the 
satellite covers approximately i of the sky, we can 
estimate that there are around 300 GRBs per year, 
for which an optical afterglow could be detected 
with timely observations. To obtain their gen- 
eral properties, we used observations published in 
the Ga mma Ray Burst Coordi nate Network Cir- 
cularfl (iBarthelmv et al.lll995l ). We chose GRBs 
detected between 2006 September and 2009 April 
that had an optical afterglow detected (and rea- 
sonably well sampled). In general, not many opti- 
cal afterglows have been detect ed in the first few 
minutes after the initial trigger ( Kann et al.ll2010l 
20111 ). Consequently, we used measured R mag- 
nitudes at approximately (t — to) =0.01 day af- 
ter the prompt GRB (tg is the time of the GRB 
trigger), at which the number of detected after- 



'^http://www.rssd.osa.int/index.php?project=GaJa& 

page=index 
^http: //gen. gsfc.nasa.gov/ 
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Fig. 2. — Number of observed GRB afterglows 
vs. their i? magnitudes shifted to 1 minute after 
the GRB trigger. The distribution is fitted with a 
normal function with mean value [x = 15.35 ± 0.09 
and standard deviation u = 1.59 ±0.09. Numbers 
are based on data from GCN circulars. 

glows is larger; if an afterglow was first detected 
at a later time, we estimated its R magnitudes 
at (i — to) =0.01 day from its light curve. Our 
sample includes 100 GRBs, with their optical af- 
terglow magnitudes at (t — to) =0.01 day spanning 
a range of values between i? = 13 mag and i? = 23 
mag. 

In general, GRB optical afterglows follow a 
power-law decay, F^{t) oc i~", with typical val- 
ues of temp oral decay index a sca ttered between 
0.4 and 1.4 (JAkerlof fc Swanll2007h . These values 
correspond to the time before a jet break, which is 
evident in some, but not all afterglows. To model 
the afterglows from early time on (where obser- 
vations are still rare), we shift the distribution of 
R magnitudes at (< — <o)=0.01 day to (t — to)=l 
minute, using the average temporal index a — 1.0. 
The time scale of 1 minute usually corresponds to 
the onset of afterglow emissiorO. We somehow ar- 
bitrarily chose the time 1 minute after the GRB 
as the start of an afterglow in order not to miss 
any possible detection. This choice does not al- 
ter our results, since, as will be shown in Section 
14.11 the probability of afterglow detection in the 
first 0.01 day after the burst is small and, there- 
fore we could have chosen any time between 1 



^ There are examples of prompt optical emission in the first 
minute after a trigger but they are scarce. 



minute and 0.01 day. We find that the distribu- 
tion of magnitudes at {t — to)=l minute is well de- 
scribed by the normal distribution with the mean 
value /i = 15.35 ± 0.09 mag and standard devia- 
tion a = 1.59 ± 0.09 mag. The distribution of R 
magnitudes, shifted to {t — to)=l minute, and cor- 
responding normal function are shown in Figure 

m 

The temporal and sky distribution of GRBs is 
random, with about 300 GRBs per year, for which 
an optical afterglow could be detected. Follow- 
ing that, in our simulation we generate 300 x 5 
GRB events distributed randomly (but uniformly) 
over the sky and randomly (but uniformly) in time 
(i.e., random to) over the five year period. Magni- 
tudes of generated GRB afterglows at [t — to)—! 
minute are distributed according to the aforemen- 
tioned normal distribution. Next we assume that 
the magnitudes are, assuming a ~ 1.0, evolving 
with time as: 

m(t~to) =m(l min) + 2.51og I ~ .° ) . (1) 

\ 1 min / 

Since the distribution of GRBs is random both in 
space and time, we do not have to specify a spe- 
cial initial position of the telescope's axes. We 
simply put the first telescope (number 1 in Figure 
[T|) at the origin of a fixed coordinate system at 
time t = and from there on compute the posi- 
tion of the telescope's axis at time t > according 
to the scanning law. After we generate a GRB, we 
first calculate the time tum after which the GRBs 
magnitude will fall below the limiting magnitude. 
We then check where the positions of both tele- 
scopes are at (i — io)=l minute, i.e., if they are 
pointing in the direction of a GRB. If they are 
not, we calculate the telescopes' axes paths over 
the sky. We stop the simulation when the time 
iiiin expires and check if there were any detections 
with the first or second telescope. After that, the 
axes of telescopes are returned to the initial posi- 
tion, another GRB is generated, and the described 
procedure is repeated. 

3.2. Orphan Afterglows 

The simulation in this part is done in a similar 
way as in the case of on-axis afterglows. Again, 
we generate a number of events, distributed ran- 
domly (but uniformly) in the sky and randomly 
(but uniformly) in time. The difference from the 



preceding section is that the OAs have not been 
conchisively observed yet, and therefore we do not 
have any observational data available on which to 
base the number of observable OAs and their light 
curve properties. We approach the problem by 
constructing a number of OAs' light curves with 
characteristics from theoretical predictions. 

Several semianalytical models have been pro- 
posed to explain the dynamics of the jet and, 
consequently, the shape of the OAs' light curve 
(|Granotl I2007T ). In our s i mulat ion, we follow 
the approach of Zou et al.l (20071 ), which in it 



self is based on lWu et al.l (|2004l ) and lGranot et al. 



(|2002l ). where an adiabatic jet with a half-opening 
angle 0j and without sideways expansion is con- 
sidered. The latter is a good approximation for 
a highly relativistic jet, as have been sh own in 
hydrodynamical simulations (lGranotll2007ll . 

With the knowledge of the jet half-opening 
angle 9j, its total collimation-corrected energy 
i^j, redshift z, and the density of the circum- 
burst medium n, we can compute the time ij 
at which an on-axis observer will see the break 
in a light curve and the flux density fj^j = 
/(tj , Ej , Z?L, z ^v,p,ee,SB) at t hat time (equations 



[2] and [3] in IZou et al.ll20071 ): Dl stands for lu- 
minosity distance, p is the power-law index of 
shock-accelerated electrons, and Eq and Eb are 
the energy equipartition factors of the electrons 
and magnetic field, respectively. The rest of the 
light curve is obtained by applying the power-law 
temporal decay with a prebreak temporal index 
ai = 1.0 (the same as in the o n-axis case) and a 
postbreak index a2 = ai + 3/4 (jZou et al.ll2007l ). 
The light curve seen by an off-axis observer is 
obtai ned by considering a point-source approxima- 



tion (Granot et al 



2002h : 

= (1 - /3)/(l - /3cos0obs) and 



(2) 



where a 

\J\ — l/r^. This approach works only in the case 
of ^obs > ^j • Ti me evolution of the Lorentz factor 
is described by (JGranot et al.ll2002l ): 



V(t) = 



-'(4) 



-3/8 



-1/2 



t < t 



t > ti 



(3) 



In Figure 13] three afterglow light curves are shown 
for the case of 0j = 0.1 and z —\. The fiux is cal- 
culated for an observed frequency v = 4.55 x lO^'' 



lc+06 

100000 

10000 

1000 

100 

10 

1 

0.1 







Z=l 


^\^^ 




ej = o.i 


- 8„b. = <i^»\ 




- 




^\ 


/ 




^^<^ 


M,i„ = 20mag 


f\^ 


■■■8„b, = 0-12...... 

..■■■■■■■C = 0-20 


'obR 


\ 





1 10 

time after t^ [days] 



Fig. 3. — Example of an optical afterglow light 
curve {i^ — 4.55 x 10^* Hz) seen by an on-axis 
observer (solid line) and the same event seen by 
an off-axis observer at an angle 9ohs = 0.12 and 
0.20 (dotted lines), to is the time at which the 
burst occurred. 



Hz. The light curve drawn with full line cor- 
responds to an on-axis observer, while the dot- 
ted lines represent the light curves seen at angles 
6»obs = 0.12 and 0.20. The lightcurves of OAs can- 
not be observed immediately after the GRB. How- 
ever, the flux slowly rises and eventually reaches 
the peak value, after which it begins to decay 
again. Depending on the parameters mentioned 
in this section, the OA's light curve can be seen 
during the period in which the OA's fiux is above 
the observer's detection limit. 

In our simulation, a light curve for a given set 
of parameters is calculated. It is then verified, 
whether it rises above the limiting flux or not. In 
the former case, the time interval tobs in which the 
OA can be observed is calculated. If during this 
time interval the position of the burst comes into 
the field of view of one of the Gaia telescopes, we 
consider that the afterglow is detected. 

Parameter values used in our simulation are: 
^j = 1 X 10^1 ergs, p = 2.2, e^ = 0.1, Eb = 0.01, 
and v — 4.55 x 10^^ Hz. We assign the same 
value of Ej to all generated bursts, since it is be- 
lieved that different observed energy values cor- 
respond to differ ent jet half-opening angle values 
(jZou et al.ll2007l ). and we do not expect that this 
will have substantial effect on our results. Other 
parameters need a more detailed discussion. 



3.2.1. Number of GRBs 

Several estim ations of the detectability of OAs 
have be en made ( Zou et al.l2007l : iTotani fc Panaitescu 
[2002^ Nakar et al.ll2002l ). Their results depend on 
jet models, initial half-opening angle distribution, 
medium surrounding the GRB source, etc. How- 
ever, we decided to make a simple estimation of 
the number of all GRBs that occur in the universe 
up to some redshift z, as described subsequently. 
During the simulation, the events with ^obs < ^j 
are sorted out and we are left only with the bursts 
for which OA light curves can be generated. 

Following [Cuetta et al.l (|2005[) the GRB forma- 
tion rate is approximated to f ollow the Rowan- 



Robi nson star formation rate ([Rowan- Robinson 
199i) 



-Rgrb — Po 



2q0.75z 
100.75 



z< 1 

z > 1, 



(4) 



wher e po ^ 33 /igg Gpc^'^ yr^^ (jGuetta et al. 
120051 ) . Considering the expansion of the universe 
(which is taken to be flat, i.e., flk = 0), the num- 
ber of bursts dN that occur between redshift z and 
z + dz in one year is: 



dNiz)^^^dz, 
^ ^ 1 + z dz 



(5) 



where the division factor (l + z) is due to time dila- 
tion. The comoving volume element dV/dz equals 



C 4:TTX^{z) 



dV _ 

dz ^ HoE{ni,z) 



X{z) 



dz 



Ho Jo E{Q,,z') 



(6) 



with E{ni,z) = ^/^i^T+YfTJlX. The number 
of GRBs per year is: 



N = 



Rgrb dV 

I + z dz 



dz. 



(7) 



In our calculation, we adopt standard cosmology 
of a flat universe with h = 0.72, £7™ — 0.3, and 
r^A =0.7 and set z,nax = 5. Thus, the calculated 
number of GRBs that occur per year is iV ^ 25000. 

3.2.2. GRB Redshift Distribution 

We generate the redshift distribution of GRB 
events according to equation ([5]) in the following 
way. The probability of findi ng a GRB at a red- 
shift z is (jCowardet al.ll2002l ): 



p{z) 



(dN/dz) 



J^"^^''{dN/dz)dz' 



(8) 



Next, a cumulative function P{z) can be calcu- 
lated as: 



P{z) 



p(z )dz , 



(9) 



giving the probability of an event occurring in the 
redshift range from z = Q io z. By calculating the 
inverse of P{z), the redshift as a function of the 
cumulative probability is obtained. Employing a 
random number generator to select values P{z), a 
GRB redshift distribution according to equation 
([5]) is generated. 

3.2.3. Jet Half-Opening Angle Distribution 

Not many GRB afterglows have been observed 
with a high enough precision to enable us the 
determination of 0j, which can be calculated as 
(iFrail et al.ll200ll ) 



6lj = 0.057 



U 



1 day 

1/8 



3/8 



-3/8 



V 1053 erg y V0.2/ VO.lcm-3/ 



(10) 



where Eiso{l) is the gamma energy output assum- 
ing a spherical burst, n^ is the efficiency in con- 
verting the energy of the ejecta into gamma rays 
and n is the number density of the circumburst 
medium. 

Using equation ([10]), we calculated angles Sj for 
43 GRB s from observation al data collected in Ta- 
bic 1 in iGao fc Pail ( 20101 ) . The obtained angle 
distribution is shown in Figure 21 In our calcu- 



lation, we set n^ 



0.2. We chose two different 



values of number density and ran the simulation 
for each number density case separately. The an- 
gle distributions have been fitted with a normal 
distribution function, and in the simulation, ran- 
domly generated half-opening angles are taken to 
be distributed according to these two functions. 
As an additional constraint we generate only an- 
gles greater than 6'j = 0.01 rad. 

Another parameter in the simulation is the ob- 
servational angle ^obs, which is distributed uni- 
formly between and -k (due to jet bimodality) . 

4. Results and discussion 

Here, we present results for both parts of the 
simulation. All results are given for a time period 
of five years, which is the expected operational 
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Fig. 4. — Distribution of calculated jet half- 
opening angles Oj for two different values of n. In 
both cases the distribution is fitted with a normal 
distribution with the mean value /i and standard 
deviation a. 

time of Gaia. All results are averaged over 1000 
and 100 simulations for on-axis and OA, respec- 
tively. 

Figure \5\ shows the number of GRB optical af- 
terglows detected by Gaia in five years depending 
on the value of Gaia^s limiting magnitude. Re- 
sults for the case of the on-axis simulation are not 
very promising. In the nominal Gaia Mnai = 20 
mag, we can expect to see less than 10 on-axis af- 
terglows in five years. Possible detection of OAs 
is more promising: with Miim — 20 mag, we ex- 
pect to detect several tens of OAs in five years. 
The latter is actually an upper limit: possible 
light extinction in the lin e of sight in the Ga laxy 



and GRB host galaxies ( Schadv et al.l 120101 ) has 



not been included in the calculation of the OAs' 
light curves. The results of the on-axis simula- 
tion, which is based on observational data, already 
include the extinction effect. Also (and this ap- 
plies to both on-axis and OAs) , the simulation was 
done for the case of the regular Johnson-Cousins R 
magnitude (Bessell 1990)), while the Gaia G mag- 
nitude is obtained with a fi lter of a much wider 
passband (j jordi et al.l 120101 ) . This results in an 
overall slightly fainter limiting magnitude. 

Hereafter, all results are given for Afii„i = 20. 

For the on-axis case, the distribution of initial 
magnitudes of optical afterglows, which were de- 
tected by (at least) one of the Gaia telescopes, is 
given in Figure [5] (left). It is a normal distribu- 





160 




140 




120 


c 

i 


100 

80 


pa 


60 

40 




20 



01-phan: n = 1 .0 cm' — • — / 

orphan: n = 0.1 cm'^ ■ P 

on-axis ▼■ / 

^r •■'' 



17 17.5 



18.5 19 



19.5 20 20.5 



21 



Fig. 5. — GRB on-axis and OA detection rate in 
five years of Gaia operation as a function of Mn^. 
In the case of OAs, the results are given for two 
different circumburst environments. 



tion with the mean value '^ 13.5 mag. This was 
expected, since the afterglows with brighter initial 
magnitudes are more likely to be detected (their 
tiim is longer). It is interesting to see the distri- 
bution of magnitudes at the time of their first de- 
tection (gray distribution in Figure [6]). Since the 
majority of detected afterglows are already quite 
faint, we cannot expect many of them to be de- 
tected by the second telescope. Indeed, the distri- 
bution of twice-detected afterglows shown in Fig- 
ure [6] in red tells us that, given the small number 
of detections with one telescope, the probability 
to detect an afterglow with both telescopes is low. 
A similar plot was made for the OA simulation 
and is shown in Figure [5] (right). The distribu- 
tion of magnitudes in the case of n = 0.1 cm~'^ 
is not given here, since it strongly resembles that 
of n — 1.0 cm~^ (but with the lower number of 
detections) . The fraction of twice-detected OAs is 
larger here. The decay in the on-axis case is con- 
stant with a = 1, while in the OA case the flux 
could actually still be rising at the time of first de- 
tection. In addition, the flux decay around Mmax 
is much shallower than a — 1, which explains the 
larger fraction of second detections in OA case. 

To check the sensibility of OA simulation re- 
sults, we look at the correlation between redshift 
z and jet half-opening angle 9j, shown in Figure 
[7] (left). As expected, most 9j values are concen- 
trated around the peak of 9j distribution (section 
13.2. 3p . Most points have values of z < 2. This can 
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Fig. 7. — Correlation between redshift z and jet half-opening angle 0j (left) and observing angle 0obs (right) 
for detected OAs. Black points represent OAs detected with one telescope, and red points represent OAs 
detected by both telescopes. 



be understood, since higher redshift value results 
in a lower observed flux. In addition, the GRB 
redshift distribution fsection F3.2.2p peaks at z = 1 
and starts to decline at higher z. Nevertheless, a 
few OAs with high redshift are detected in the sim- 
ulation. Since all of these OAs have relatively low 
0j , which neutralizes the negative redshift effect on 
the flux value, their detection is understandable. 
The majority of second detections have low red- 
shift values and no preferable By This is mostly 
because the majority of detections with the first 
telescope have low redshift values. Because of the 
relatively flat light curve around the peak value 
(Figure [3]) , once the OA has been detected with 



the first telescope, redshift and 0j are not the cru- 
cial parameters for the second detection. In order 
to detect OA with the second telescope, the time 
t — io of the first detection, i.e., before or after the 
peak flux value, is important. 

Figure [7] (right) also shows the correlation be- 
tween redshift z and observing angle ^obs (which 
for the OAs has to be larger than Q-^ of a particu- 
lar burst). If a burst has a low redshift, 0obs can 
be quite large, as is evident in the long tail. On 
the other hand, a very low observing angle helps in 
the detection of an OA from a distant object, since 
the difference between 6'j and 6'obs is small and the 



burst is 'almost' on-axis. Detections with the sec- 
ond telescope do not appear to have a preferable 
position in the graph. 

4.1. Identification of GRB Optical After- 
glows 

On-axis afterglows are initially rapidly fading 
and they will be observed for 4.4 s by each of 
nine transiting CCD detectors. In principle, the 
change of magnitude in that time, i.e., in observa- 
tions by two successive CCDs or CCD field transit, 
could be observable. Therefore, we calculated the 
change of magnitude in 4.4 s for all detected on- 
axis afterglows and also for observations lasting as 
long as a particular afterglow is in the telescope's 
field of view. The distribution of the change of 
magnitude for both cases is shown in Figure |8] 
(left). We calculated the changes only for the 
first afterglow detection, since the probability of 
detection with the second telescope is small. In 
practice, the change of magnitude could be de- 
tected only if it was larger than the photometric 
error. The expected photometric errors, averaged 
across the sky, in the case of Mn^ = 20 mag will 
be 10 mmag at bes t (though it will be s maller 
for brighter events) (jPerrvman et al.l 12001 ). The 
chances of observing the change in magnitudes will 
thus be limited. The change in magnitude of OAs 
will be even harder to detect, since the flux rises 
and decays more slowly. Also, OA could actu- 
ally become brighter through the observation if the 
detection happened prior to the peak light-curve 
value (Figure m right). 

Looking at the time of detection relative to the 
initial GRB (Figure O, most of on-axis afterglows 
are expected to be detected at ^ 0.1 day after 
the GRB. The time of detection is considerably 
larger for OAs. Since they are still bright a few 
days after a GRB, in addition to a considerable 
probability of prepeak detection, they could also 
be observed with Goia's second telescope and, if 
identified quickly, with ground-based telescopes. 

Since the change in afterglow magnitude while 
being scanned by Gaia is not expected to be 
large enough to enable reliable way of identifica- 
tion of a transient source as a GRB optical af- 
terglow, we consider other possibilities for GRB 
afterglow identification from one short detection. 
Since Gaia will have a photometric instrument on- 
board, there is the possibility to identify after- 
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Fig. 8. — Left: Distribution of the change in mag- 
nitude in 4.4 s (corresponding to an observation 
with two successive CCDs) (black) and the change 
of magnitude in the time in which an afterglow is 
in the entire field of view (red) for the on-axis case. 
Right: Distribution of the change in magnitude in 
the time in which an orphan afterglow is in the 
field of view. Black distribution corresponds to 
detection after the peak in a light curve, and red 
distribution corresponds to detection before the 
peak. 
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Fig. 9. — Distribution of first detection times 
tdet relative to the time to of the initial GRB for 
the case of on-axis (red) and orphan afterglows 
(black). 



glow from its spectral energy distribution (SED). 
It has been observed that the SED of afterglows 
follows a power- law (F^ oc i'~f) modulated by ex - 
tinction ( Schadv et al.l l2010t iGreiner et al.l 12011 ) . 
We show the expected SED shape in UV-to-NIR 



frequency range in Figure I10[ where an average 
spectral index value /3 = 0.6 has been assumed 
(JGreiner et al.ll201l[ ). As shown in many statisti- 
cal analyses, in the majority of cases the extinction 
in GRB host galaxies is best described with a so- 
called SMC extinction profile (Pei 1991). Thus, we 
added SMC-type extinction to the SED (dashed 
lines in Figure [101), assuming Ay = 0.3, which 
is the av erage determined extinction value in th e 
V band (jSchadv et al.lboiol: [Greiner et al.ll201ll) . 
The extinction in the host-galaxy rest frame has 
been calculated for an average observed GRB red- 
shift z — 2.2 Ifl and for redshift z — 0.7, which is 
the average redshift of detected OAs in our simula- 
tion. There is also a contribution to extinction by 
our Galaxy. We left the contribution of Galactic 
extin ction to vary between < Av,ga < 0.3 and, 
using ICardelli et al.l (|l989l) , added Galactic con- 
tribution to the SED. In Figure [TUl shaded areas 
mark the possible range of SED curves depend- 
ing on the Av,GA value. Vertical lines correspond 
to the central wavelengths of the Gai a photomet- 
ric in strument's three wide passbands (jjordi et al.l 
20101 ). Since Gaia photometry will be obtained by 



means of two low-resolution spectra (so-called red 
photometer (RP) and blue photometer (BP) spec- 
tra with central frequencies shown in Figure [10] as 
Grp and Gbp vertical lines), these spectra could 
provide us with enough information to determine 
the SED shape. In the case of bright sources, in- 
tegrated flux from the radial velocity instrument 
(Grvs) can also be obtained. As shown in Figure 
[TOl the change in flux density in the observational 
spectral range (marked as AF) is expected to be 
large enough to be measurable. If the SED of a 
detected transient source resembles the one shown 
in Figure [TUl the source could be considered to be 
a GRB afterglow. 

Another possibility to identify an on-axis af- 
terglow is to cross-check the position and time of 
Gaia detection with GRB triggers by gamma-ray 
satellites (e.g., Swift, Fermi, Integral). Ground- 
based follow-up observations would also be helpful 
in afterglow recognition, but their usefulness will 
critically depend on the delay between the Gaia 
detection and issuing of the alert to the commu- 
nity (which is expected to be about 24 hr) and the 
limiting magnitude of the follow-up campaign. 
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Fig. 10. — Expected spectral energy distribution 
of detected GRB optical afterglows without ex- 
tinction (solid line), with added host-galaxy ex- 
tinction (dashed lines) and with additional Galac- 
tic extinction, where different values of extinction 
are considered (filled areas). Host-galaxy extinc- 
tion was calculated at redshift z = 0.7 (red), cor- 
responding to the average redshift of detected OAs 
in our simulation, and z = 2.2 (black), correspond- 
ing to the average observed redshift of GRBs, de- 
tected with the Swift satellite. Central frequencies 
of Gaia^s photometric bands are marked with ver- 
tical lines. See the text for more details. 



5. Conclusions 

We have shown that the number of expected 
GRB optical afterglow detections with Gaia satel- 
lite during its five years of expected operation 
is of the order of a few tens. Most of the ex- 
pected detections are of the orphan afterglows, 
while the possibility of an on-axis afterglow detec- 
tion is small (less than 10). The results are depen- 
dent on a set of various parameters, especially in 
the case of OAs. Since we made some assumptions 
and simplifications in the course of our simulation, 
we would like to point out that these results have 
to be interpreted as an order-of-magnitude estima- 
tions. In addition to the small number of expected 
detections, poor time sampling of the afterglow 
light curve and problem of afterglow identification 
make the study of GRB optical afterglows with 
Gaia difficult. On the other hand, since OAs have 
not been conclusively detected yet, a possible de- 
tection would provide valuable new data. 



^ http://heasarc.gsfc.nasa.gOv/docs/5tui/f/archive/grb_table/ 
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A. Gaia's scanning law 

To describe Gaia^s scanning law we first set aside its rot ation around the Sun . Angles 7, ftt and ojt (at 



an arbitrary time t) represent the well-known Euler angles (jHand fc Finchl 119981 ). Coordinates of the first 
and second telescope axes positions in the fixed coordinate system (5') with its origin at the satellite are 
described with ri^2, while coordinates in the system of the scanning satellite (S") are described with r'i^2- 
Transformation between those two systems is given by matrix U 

Via = Uv\a, (Al) 

which is obtained by multiplying three rotational matrices, each one describing a rotation with one of the 
Euler angles: 

/ cos(a;t) cos(f2i) — cos(7) sin(f2t) sin(a;t) — sin((.ijt) cos(f2t) — 003(7) sin(f2t) cos([.i;t) sin(7) sin(r2t) 
U = cos(ajt) sin(f!t) + 005(7) cos(f2t) sin(tjt) — sui{uit) sin(f2t) + 005(7) cos(f2t) cos{uit) — 8111(7) cos(04) 
y sin(7) sin((.jt) sin(7) cos(tji) 003(7) 

The first telescope axis is placed to be the x' axis, i.e., r'l — (1,0,0). Considering the axes' separation 
angle /3, the direction of the second telescope axis is described by r'2 = (cos /?, — sin/3, 0). Now we introduce 
another (final) fixed coordinate system S*, initially coinciding with S. The rotation around the Sun (for an 
angle depending on the time of the year) is obtained by rotating S around the x axis (Figure [1]) . 
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